Neutrino Mass Hierarchy and Stepwise Spectral Swapping of Supernova Neutrino 

Flavors 



Huaiyu Duan, 1 George M. Fuller, 1 J. Carlson, 2 and Yong-Zhong Qian 3 

'Department of Physics, University of California, San Diego, La Jolla, CA 92093-0319 

2 Theoretical Division, Los Alamos National Laboratory, Los Alamos, NM 87545 
3 School of Physics and Astronomy, University of Minnesota, Minneapolis, MN 55455 

(Dated: February 1, 2008) 

We examine a phenomenon recently predicted by numerical simulations of supernova neutrino 
flavor evolution: the swapping of supernova u e and z^, T energy spectra below (above) energy Ec 
for the normal (inverted) neutrino mass hierarchy. We present the results of large-scale numerical 
calculations which show that in the normal neutrino mass hierarchy case, Ec decreases as the 
assumed effective 2x2 vacuum v £ ^ ^ MlT mixing angle (~ #13) is decreased. In contrast, these 
calculations indicate that Ec is essentially independent of the vacuum mixing angle in the inverted 
neutrino mass hierarchy case. With a good neutrino signal from a future Galactic supernova, the 
above results could be used to determine the neutrino mass hierarchy even if #13 is too small to be 
measured by terrestrial neutrino oscillation experiments. 

PACS numbers: 14.60.Pq, 97.60.Bw 



In this letter we point out how two grand themes 
in contemporary science, physics beyond the Standard 
Model of elementary particles and the physics of stars 
undergoing gravitational collapse, overlap in a way that 
could allow unique insight into the nature of neutrinos. 
Recent experiments have established that neutrinos have 
non- vanishing rest masses and that the flavor states u e , 
Vfj, and v T for these particles are mixtures of the vacuum 
mass eigenstates v±, v-i and 1^3 (see, e.g., Refs. [HQ for 
recent reviews). However, key issues remain unresolved. 
Among these is the nature of the neutrino mass hier- 
archy: the sign of the mass-squared difference Sm^ = 
to 2 — m\ ~ ±(Sm 2 tm remains unknown. Here Sm1 tm is 
the neutrino mass-squared difference associated with at- 
mospheric neutrino oscillations, and the plus (minus) sign 
corresponds to the normal (inverted) neutrino mass hier- 
archy. Conventional laboratory experimental resolution 
of the mass hierarchy issue is problematic, in part because 
neutrino rest masses are tiny and because 613, the mix- 
ing angle relating v e to ^3 , is small. One possible way to 
probe the neutrino mass hierarchy is to analyze neutrino 
signals from Galactic supernovae. Supernova neutrinos 
can experience significant flavor transformation through 
the Mikheyev-Smirnov-Wolfenstein (MSW) effect d [| 
as they stream out from the surface of a proto-neutron 
star (with very high matter density) into the vacuum. In 
addition, it has been pointed out that neutrino-neutrino 
forward scattering can provide an additional source for 
neutrino refractive indices 0, S S Hi • This neutrino self- 
coupling is especially important in the supernova envi- 
ronment because neutrino fluxes are large. Therefore, 
previous studies of supernova neutrino oscillations based 
on the pure MSW effect % M, EI El, El EI El may 
not apply in some scenarios. Following suggestions in 



Refs. jig, [17|, [18( , in this letter we present and analyze 



new large-scale numerical calculations of supernova neu- 



trino flavor evolution that suggest a novel method to de- 
termine the neutrino mass hierarchy. This method is 
independent of absolute neutrino masses and can work 
even for tiny #13. 

Our method is based on a stunning feature revealed by 
recent numerical simulations of supernova neutrino flavor 
evolution 16|, HJ: (1) For the normal neutrino mass hi- 
erarchy case, v e and f M/r swap their energy spectra at 
energies below a transition energy Ec, but retain their 
original spectra at higher energies; (2) For the inverted 
neutrino mass hierarchy case, the situation is exactly the 
opposite. This phenomenon is known as "stepwise spec- 
tral swapping" (lij or "spectral split" [l8j |. 

The stepwise swapping of the v e and v^ T energy spec- 
tra has its origin in nonlinear neutrino self-coupling. As- 
suming coherent neutrino propagation and the efficacy of 
the mean field approach 2^, 21 1, for 2 x 2 flavor evolu- 
tion, it is possible to define the neutrino flavor isospin 
(NFIS) as [22| 



and s p = (o-y-ilivf — (o-y^v) (1) 



for a neutrino (with flavor wavefunction Vv) and an an- 
tineutrino (ipa), respectively. Flavor evolution for neu- 
trino or antineutrino mode i is described by precession 
of a corresponding NFIS Sj around an effective field: 
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Here el „ _ are the flavor-basis unit-vectors in flavor 



space, Hy = — sin 29 v e l x + cos 26 v e l z generates vacuum 
mixing for a nonvanishing effective mixing angle 9 V , 
oji = ±5m 2 /2Ei is the vacuum precession angular veloc- 
ity around Hy for a NFIS corresponding to a neutrino 
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(plus sign) or antineutrino (minus sign) with energy Ei, 
Gf is the Fermi constant, n e is the net electron number 
density, is the angle between the propagation direc- 
tions of neutrinos in modes i and j, and rij is the number 
density of neutrinos in mode j. Because flavor transfor- 



mation in the 



and 



channels is the 



most important in supernovae (e.g., for shock reheating 
HHl and nucleosynthesis ^QQlil), and because 



Sm 2 



will give flavor transformation deeper in the su- 



pernova envelope than will Sitiq , the mass-squared differ- 
ence associated with solar neutrino oscillations, we take 
Sm 2 = ±3 x I(T 3 eV 2 ~ ±Sml tm and V ~ 6 13 < 1. For 
this 2x2 mixing we use v T * to designate the relevant 
linear combination of and v T [28] . 

In a stepwise-swapping scenario, the probability P vv 
for neutrinos to remain in their initial flavor states is 



P w M~i[l-sgn(u,-a,° r )], (3) 

where sgn(£) = £/|£| is the sign of £, and Wp r = 5m 2 /2Ec 
specifies the transition energy Eg- Eq. ([3]) is slightly 
different from Eq. (57b) in Ref. [29(, with different con- 
ventions for V and Sm 2 . This stepwise spectral swap- 
ping feature has been demonstrated in two different ap- 
proaches 0, 0, 31 j : "multi-angle" simulations, 



where flavor evolution on independently-followed neu- 
trino trajectories is self-consistently coupled, and "single- 
angle" simulations, where the evolution history of ra- 
dially propagating neutrinos is assumed to apply to all 
trajectories. Single-angle calculations capture the qual- 
itative features of stepwise spectral swapping, and they 
suggest the following generic explanation for this phe- 
nomenon. 

Because neutrinos and antineutrinos are in flavor 
eigenstates when they leave the neutrino sphere, they 
naturally form a "bipolar system" in which the corre- 
sponding NFIS's form two oppositely oriented groups 
|2a |. (Note that v e /v T * and v e jv r * correspond to NFIS's 
in the directions of +e l z and — e { z , respectively.) This 
bipolar system behaves like a gyroscopic pendulum in 
flavor space (32] |. which can have both nutation and pre- 
cession modes. These neutrinos and antineutrinos ini- 
tially follow a quasi-static, MSW-like solution near the 
neutrino sphere before being driven away from this so- 
lution by the collective nutation of the gyroscopic pen- 
dulum [291] . Subsequently, the gyroscopic pendulum can 
execute regular precession around Hy, corresponding to 



the collective precession of NFIS's in flavor space [17], 
This precession, although not perfectly regular, is indeed 
found in both single-angle and multi-angle simulations 

Hill]. If NFIS's stay m the regular collective preces- 
sion mode, a stepwise spectral swapping given by Eq. ([3]) 
will occur when the neutrino fluxes decrease toward 

3 12, HI- In this case, aA in Eq. ([3]) is just the pre- 



cession angular velocity for vanishing neutrino fluxes. 
Strictly speaking, the regular precession mode obtains 



for n e = where the "lepton number" 

/>oo 

/ [f Vl (E)-f„ 3 (E)-fc(E)+f P3 (E)}dE (4) 
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is conserved j!!]. Here f Vl / Pl \(E) and /i/ 3 (p 3 )(-E0 are 
the distribution functions specifying the populations of 
the corresponding neutrino (antineutrino) vacuum mass 
eigenstates within energy interval dE. These are normal- 
ized to the total (summing over all states) neutrino num- 
ber density n* ot and in general evolve with time. The 
conservation of C is exact for n e — 0, and holds even 
when neutrino number densities change [171 ] . Because 
the presence of the matter field does not change the col- 
lective precession qualitatively 1171 12911 . the conservation 
of C can be used to compute Eq [18| . 

In the inverted neutrino mass hierarchy case, flavor 
transformation is suppressed when n e and neutrino fluxes 
are high. As v e are dominant in supernovae, the bipolar 
system of neutrinos and antineutrinos resembles a gy- 
roscopic pendulum near its highest point (displacement 
angle equal to 7r) in flavor space. For a simple bipolar 
system initially consisting of mono-energetic v e and v e , 
the analogy is exact, with the initial displacement angle 
being 7r — 29 v . When the total neutrino flux decreases 
below some critical value, the flavor pendulum evolves 
away from its maximum displacement 17(, |32j. Its nu- 
tation then pushes neutrinos and antineutrinos into the 
collective precession mode. The presence of a matter 
field does not change this nutation qualitatively, but ef- 
fectively reduces the mixing angle [i^, H2]. Therefore, 
C is essentially unchanged before the precession mode 
begins. As a result, when n e and neutrino fluxes be- 
come small, fuy{ui){E) and f V3 (v 3 )(E) are related simply 
through Eq. (|3|) to the initial neutrino energy spectra at 
the neutrino sphere. Specifically, for 9 V <C 1 we have 
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where, e.g., f Va {E) and n V(s (E) are the initial spectrum 
(normalized to n* ot ) and number density, respectively, of 
v e at the neutrino sphere. The transition energy Eq can 
then be found from Eq. ([5]). 

The conservation of C can also be used to find Eq for 
the normal mass hierarchy case. However, in this case, 
C cannot be related simply to the initial neutrino spec- 
tra at the neutrino sphere through Eq. ([3]). This is be- 
cause there is a resonance in the quasi-static MSW-like 
solution initially followed by neutrinos and antineutri- 
nos. For example, in the large neutrino luminosity limit, 
all neutrinos and antineutrinos are synchronized [331 ] and 
experience simultaneously an MSW-like resonance near 
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FIG. 1: (Color online) The neutrino survival probability P vv as a function of neutrino emission angle $o (relative to the normal 
at the emission point on the neutrino sphere) and energy E. The left panel is calculated for a normal neutrino mass hierarchy 
with 8 V — 0.01 and the right panel is for an inverted neutrino mass hierarchy with 9 V = 10 -9 . These results are taken at radius 
r = 250 km. Except for 8 V , all parameters are the same as those for Fig. 3 of Ref. [l9j . 



the radius where a single v e with a representative en- 
ergy i^svnc would encounter a conventional MSW reso- 
nance [26|. During this MSW-like evolution, C is de- 
creased. Using the initial spectra for supernova neutrinos 
and Eqs. ([2]) and (01, it can be shown that the less C is 
reduced, the smaller Eq becomes. If 9 V is tiny and/or the 
neutrino luminosities are not large enough, the MSW-like 
conversion will be non-adiabatic, and C will change very 
little. In this case, Ec — > and the stepwise nature of 
the swapping of v e and i^ 1jT spectra becomes unobserv- 
able. On the other hand, given sufficiently large values 
of 9 V and/or neutrino luminosities, MSW-like flavor con- 
version can be adiabatic and efficient for a large range 
of neutrino energies. This can engender nearly complete 
swapping of the entire neutrino and antineutrino spectra 
and, consequently, L can retain its magnitude but reverse 
its sign. In this case, Ec for the subsequent stepwise 
spectral swapping is roughly the same as in the inverted 
neutrino mass hierarchy case. 

Though the transition energy Eq is sensitive to 9 V in 
the normal mass hierarchy case, it appears to be essen- 
tially independent of 9 V in the inverted mass hierarchy 
case. We have carried out multi-angle simulations under 
the same conditions as those discussed in Refs. 
except with smaller 8 V . The probability P l , cU<l (E,'do) at 
radius r = 250 km is plotted in Fig. [T] for both a normal 
mass hierarchy case (9 V = 0.01, left panel) and an in- 
verted mass hierarchy case (9 V — 10" 9 , right panel). Here 
■9q is the angle between the propagation direction of the 
neutrino and the normal at its emission position on the 
neutrino sphere (see Fig. 1 of Ref. [Hj|). A comparison of 
the results shown in Fig. [1] with those for 6 V = 0.1 shown 
in Fig. 3 of Ref. [19j is revealing. In the normal neutrino 
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FIG. 2: Single-angle simulation results for the dependence of 
rx on 9 V in the inverted neutrino mass hierarchy case. Here 
9 V is the effective 2x2 vacuum mixing angle (~ 613), and rx 
is the radius where the energy- averaged value of P VcV<1 drops 
below 0.9. Except for 9 V , all parameters are the same as those 
for Fig. 2 of Ref. \v^. 

mass hierarchy case, Eq decreases from ~ 10 MeV to ~ 3 
MeV as 9 V is reduced from 0.1 to 0.01. However, for the 
inverted neutrino mass hierarchy case, Ec is essentially 
unchanged as 9 V is decreased by 8 orders of magnitude. 
The value E c ^ 8.4 MeV calculated from Eq. |JSj! for this 
case agrees very well with the numerical results. We note 
that Ec has a slight dependence on neutrino trajectory 
(cos#o) in the multi-angle simulations for 9 V -c 0.1. 

The apparent insensitivity of Ec to 9 V in the inverted 
neutrino mass hierarchy case requires discussion. We 
note that the neutrino system transitions from the MSW- 
like evolution to the collective precession mode through 
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nutation. If the system does not develop significant nu- 
tation while it is in the collective flavor transformation 
regime, it will not enter the collective precession mode, 
and therefore, stepwise spectral swapping will not oc- 
cur. For a uniform and isotropic gas of mono-energetic 
neutrinos initially in pure v e and v e states, the nutation 
timescale is T nut ~ — In 8 V in the inverted mass hierar- 
chy case [23. Estimating the nutation timescale for re- 
alistic supernova neutrino systems is problematic, partly 
because this quantity depends on time- varying electron 
and neutrino number densities 22, [32|, [34[. In Fig. [2] 
we plot as a function of V the radius rx (as defined in 
Ref. [Til]) where the energy- averaged value of -P„ cl , e drops 
below 0.9 and significant nutation develops in our single- 
angle simulations for the inverted mass hierarchy. Our 
single-angle calculations suggest that the onset of signif- 
icant nutation is nearly independent of 6 V as it is de- 
creased from 0.1 to ~ 10~ 20 . As 9 V is decreased further, 
rx begins to increase. We expect that for sufficiently 
small 6 V , the onset of significant nutation is pushed to so 
large a radius that the corresponding neutrino number 
density becomes too low to generate any collective flavor 
evolution or stepwise spectral swapping. 

Just after the bounce of the supernova core, when 
the supernova shock breaks through the neutrino sphere, 
there is a brief intense burst of neutrinos which are emit- 
ted predominantly as v^s. Lacking in D^s [351] . this burst 
is not likely to be affected by stepwise spectral swapping 
at the Sm1 tm scale (but see Ref. 36]). Later, both neu- 
trinos and antineutrinos are emitted and they form a 
bipolar system. If #13 is not too small, for the normal 
mass hierarchy case, stepwise spectral swapping at the 
Sm^ tm scale can be observable at late times when matter 
has been sufficiently condensed toward the proto-neutron 
star. For the inverted mass hierarchy case, however, step- 
wise spectral swapping occurs even at early times [22, 31] 
for essentially any nonvanishing 6*13. In light of the in- 
sensitivity of Eq to #13 in the latter case, supernova neu- 
trino signals can offer a unique probe of the neutrino 
mass hierarchy even for #13 too small to be measured by 
conventional neutrino oscillation experiments. 
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